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Abstract 

Lattice QCD with Wilson quarks and a chirally twisted mass term (tmQCD) 



D 
> 

has been introduced in refs. |ljj2j. We here apply Symanzik's improvement 
programme to this theory and list the counterterms which arise at first order 
in the lattice spacing a. Based on the generalised transfer matrix, we define 
the tmQCD Schrodinger functional and use it to derive renormalized on-shell 
correlation functions. By studying their continuum approach in perturbation 
theory we then determine the new 0(a) counterterms of the action and of a 
few quark bilinear operators to one-loop order. 



1 Introduction 



In ref. 0] twisted mass lattice QCD (tmQCD) has been introduced as a solu- 
tion to the problem of unphysical fermion zero modes which plague standard 
lattice QCD with quarks of the Wilson type. We will assume that the reader 
is familiar with the motivation of this approach, and refer to Q for an intro- 
duction. The main topic of the present paper is the application of Symanizik's 
improvement programme to tmQCD. We introduce the set-up in the simplest 
case of two mass-degenerate quarks, and study the improved action and the 
improved composite fields which appear in the PCAC and PCVC relations. 

Our strategy follows closely refs. in section 2 we go through the 

structure of the 0(a) improved theory. We then define the Schrodinger func- 
tional for tmQCD, and use it to derive suitable on-shell correlation functions 
(sect. 3). The perturbation expansion is then carried out along the lines of 
ref. ||], and the new 0(a) improvement coefficients are obtained at the tree- 
level in sect. 4 and to one-loop order in sect. 5. A few details have been 
delegated to appendices. Appendix A describes how the twisted mass term 
can be incorporated in Liischer's construction of the transfer matrix ||, and 
Appendix B contains the analytic expressions for the coefficients used in the 
analysis of the one-loop calculation. 

2 Renormalized and O(a) improved tmQCD 

The renormalization procedure for twisted mass lattice QCD with Wilson 
quarks has already been discussed in ref. ||]. Here we apply Symanzik's im- 
provement programme to first order in the lattice spacing a. The procedure 
is standard and the details of its application to lattice QCD with Nf mass 
degenerate Wilson quarks can be found in ref. Q. 

Our starting point is the unimproved tmQCD lattice action for a doublet 
of mass degenerate quarks, 

S[U,$,1}] = S g [U\ + Sf[U,$,iI>], (2.1) 

with the standard Wilson gauge action and the fermionic part 

S F [U,ij,^] = a 4 ^(x) (D + m + i W T 3 )^(x). (2.2) 

X 

The massless Wilson-Dirac operator is given by 

D = 5 E{( V m + V t) ^ " «V*V M }> ( 2 -3) 
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where the forward and backward covariant lattice derivatives in direction \jl 
are denoted by V M and respectively. As tmQCD with vanishing twisted 
mass parameter /x q reduces to standard lattice QCD we expect that improve- 
ment is achieved by using the standard 0(a) improved theory and adding the 
appropriate 0(a) counterterms which are proportional to (powers of) /x q , and 
which are allowed by the lattice symmetries. The procedure hence consists in a 
straightforward extension of ref. ||] , and we take over notation and conventions 
from this reference without further notice. 

2.1 Renormalized O(a) improved parameters 

Following ref. Q we assume that a mass-independent renormalization scheme 
has been chosen, and we take the same steps as done there for standard lattice 
QCD. At fi q = the Sheikholeslami-Wohlert term suffices to improve the 
action, up to a rescaling of the bare parameters by terms proportional to the 
subtracted bare mass m q = tuq — m c ||. At non- vanishing /i q we find that 
improved bare parameters are of the form 

9o = 9o(l + b g am (l ), (2.4) 
m q = m q (l + o m am q ) + 6 m a^ q , (2.5) 
Jl q = /x q (l + 6 M am q ), (2.6) 

i.e. there exist two new counterterms with coefficients 6^ and b m . The renor- 
malized 0(a) improved mass and coupling constant are then proportional to 
these parameters, viz. 

gi = alz g {gi,a^), (2.7) 

m R = m q Z m (^,a/i), (2.8) 
^ R = Jl q Z^(gQ,afx). (2.9) 

The ratio of the appropriately renormalized mass parameters determines the 
angle a which is involved in the physical interpretation of the theory ]|]. We 
will discuss below the general 0(a) improved definition of a. Here we note 
that the case of particular interest, a = vr/2, corresponds to m R = 0, which 
implies m q = 0(a) 0. In this case all the usual 6-coefficients multiply terms 
of 0(a 2 ) and are thus negligible in the spirit of 0(a) improvement. One then 
remains with a single coefficient 6 m , which compares favorably to the situation 
in standard lattice QCD where two coefficients, b m and 6 g , are required. 
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2.2 Renormalized O(a) improved composite fields 

We assume that composite fields are renormalized in a mass-independent 
scheme, and such that the tmQCD Ward identities are respected |2j. At- 
tention will be restricted to the quark bilinear operators which appear in the 
PCAC and PCVC relations. Moreover, we only consider the first two flavour 
components, and thus avoid the renormalization of power divergent opera- 
tors such as the iso-singlet scalar density 0. As explained in ref. ||], the 
third flavour component of the PCAC and PCVC relations can be inferred in 
the continuum limit, by assuming the restoration of the physical isospin sym- 
metry. The 0(a) improved currents and pseudo-scalar density with indices 
a, b € {1, 2} are then parameterised as follows, 

(A R )J = Z A (l + b A am q )\A^ + c A ad^P a + a^b A e 3ab V^, (2.10) 



(2.11) 



(Vr)1 = Z V (1 + byam^) |t£ + cyad u T^ + a/x q 6 V e 6ab 
(P R ) a = Z P (1 + b P am q )P a . (2.12) 

Here we have chosen the bare operators which are local on the lattice, with 
the conventions of ref. ||. While this is the simplest choice, we also recall the 
definition of the point-split vector current, 

TO = \{${x){^-l)^U{x,ii)i>{x + afL) 

+ $(x + oA)(7a. + l)f^,|i)"V(4 (2.13) 

which is obtained through a vector variation of the action. This current is 
protected against renormalization, and the PCVC relation 

d;V«(x) = -2^ ah P h {x), (2.14) 

is an exact lattice identity, with the local pseudo-scalar density and the back- 
ward derivative d* in //-direction Q. This implies the identity Z^Zp = 1 in 
any renormalization scheme which respects the PCVC relation. 

2.3 An alternative definition of the improved vector current 

An alternative renormalized improved current can be obtained from the point- 
split current ( 2.13| ). For this it is convenient to start from the symmetrized 



version 



W = \ ( W + v > - a m > ( 2 - 15 ) 
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which behaves under space-time reflections in the same way as the local vector 
current. The counterterm structure then is the same as in eq. (2.11), i.e. one 
finds 

(V50£ = Z v(l + M™q) [K + OvaB„T^ + bya^e^A^] , (2.16) 

where we have again restricted the indices a, b to the first two components. 
One may now easily show that 

Zy = 1, b v = 0. (2.17) 

To see this we first note that at /i q = the vector charge of this current is 
given by 

Q y {t) = \Zy(\ + byam q ) [Q v (t) + Q v (t - a)} , (2.18) 

with 

Qv(*o) = a 3 ^y o a (x). (2.19) 

X 

At /U q = 0, correlation functions of the charge are xo-independentQ, and the 
0(a) improved charge algebra for and the exact charge algebra for Qy 
together imply that the whole renormalization factor in eq. (pnsl) must be 
unity. As this holds independently of m q , one arrives at the conclusion ( 2.17] ). 

A further relation is obtained by noting that the PCVC relation between 
the renormalized 0(a) improved fields, 

S M (Vk)« = -2^ R e 3ab (P R ) b , (2.20) 

with the symmetric derivative = \(d^ + d*) must hold up to 0(a 2 ) correc- 
tions. Then, using the identity 

W(x) = o; (yfc) + \a 2 d;dX{x)) , (2.21) 

one obtains the relation 

Z P Z m Z^by = -(b li + b P ). (2.22) 

The scale-independent combination of renormalization constants multiplying 
by is determined by axial Ward identities |§, so that eq. ( [2.22 ) can be consid- 
ered a relation between improvement coefficients. 

1 i.e. as long as the time ordering of the space-time arguments in the given correlation 
function remains unchanged. 
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2.4 O(a) improved definition of the angle a 

The physical interpretation of the correlation functions in tmQCD depends on 
the angle a, which is defined through 

tana = — . (2.23) 
m R 

In this equation fi R and m R are the 0(a) improved renormalized mass parame- 
ters which appear in the PCAC and PCVC relations [||. Up to terms of 0(a 2 ) 
we then find 

/% _ jU q [l + (bg - b m )am q } _ fi q [l + (6 M + bp - b A )am q ] ^ ^ 

m R Z P Z m [m q + 6 m a/x 2 ] Z A [m + ^a// 2 ^ 1 ] 

Here, to denotes a bare mass which is obtained from some matrix element of 
the PCAC relation involving the unrenormalized axial current + c^d^P 1 
and the local density P . Given m, the critical mass m c , and the finite renor- 
malization constants Z\, Zy and ZpZ m , the determination of the 0(a) im- 
proved angle requires the knowledge of two (combinations of) improvement 
coefficients, which may be chosen to be 6 M — 6 m and 6 m , or b^ + bp — 6a and 6a- 
A special case is again a = vr/2, which is obtained for vanishing denominators 
in eq. ( |2.24| ). For this it is sufficient to know either 6a 

or 6 m , and the finite 
renormalization constants Za or ZpZ m are then not needed. 



2.5 Redundancy of improvement coefficients 

Having introduced all 0(a) counterterms allowed by the lattice symmetries, 
it is guaranteed that there exists a choice for the improvement coefficients 
such that 0(a) lattice artefacts in on-shell correlation functions are completely 
eliminated. We now want to show that there is in fact a redundancy in the 
set of the new counterterms introduced so far, i.e. the counterterms are not 
unambiguously determined by the requirement of on-shell improvement alone. 
To see this we consider the renormalized 2-point functions 

G A (x-y) = ((Ar^xXPr) 1 ^)}, (2.25) 
G v (x-y) = ((Ur^xXPr) 1 ^)), (2.26) 

of the renormalized 0(a) improved fields defined in subsect. 2.2. We assume 
that a quark mass independent renormalization scheme has been chosen, and 
with the proper choice for the improvement coefficients one finds, 

G x (x) = lim G x (x) + 0(a 2 ), X = A, V, (2.27) 
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provided that x is kept non-zero in physical units. If the new improvement 
coefficients b m , 6 M , 6 a and 6y were all necessary any change of 0(1) in these 
coefficients would introduce uncancelled 0(a) artefacts in eq. ( gjg ). Varying 



the coefficients b m —* b m + A6 m , 6 M — > 6 M + A6 M and 6a — ► 6a + A6a in the 
correlation function G A (x), we find that the correlation function itself changes 
according to 

r ~ d 
AGa(i) = -a^ R Z P Ab m Z P Z m /j, R - — G A (x) 



dm 



R 



d 

+ A6 At (Z P Z m )" 1 m R ^ — G A (x) 
-AMaVG v (4 (2.28) 

where terms of 0(a 2 ) have been neglected. In the derivation of this equation 
one has to be careful to correctly take into account the counterterms propor- 
tional to 6^ and 6 m . First of all we notice that changing an 0(a) counterterm 
can only induce changes of 0(a) in the correlation function. For instance, the 
equation 



G A (x) =G A (x) + Ab ll ^-G A (x) + 0(a 2 ), (2.29) 



holds even for finite changes A6^. Second, when taking the continuum limit 
the bare mass parameters become functions of the improvement coefficients 
such that the renormalized 0(a) improved masses are fixed. For instance one 
has 

H q = Z P fi R (l - 6 jU Z~ 1 am R ) + 0(a 2 ), (2.30) 
and a straightforward application of the chain rule leads to 

where we have used eq. ( |2.30D and neglected terms of 0(a 2 ). Proceeding in the 
same way for the variation with respect to 6 m , and changing to renormalized 
parameters /x q = Zp/j, R + 0(a), m q = Z~^m R + 0(a) eventually leads to 



eq. O). 

At this point we recall eq. (3.13) of ref. ||], which expresses the re- 
parameterization invariance with respect to changes of the angle a. In terms 
of the above correlation functions one finds, up to cutoff effects, 

^-Ga(x) ee LI - MR^y G A (x) = -Gy(x). (2.32) 
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As a consequence not all the terms in eq. fl2.28|) are independent, and the 
requirement that AGa(x) be of order a 2 entails only two conditions, 

Ab m + Ab^ZpZ^- 2 = 0, (2.33) 
Ab m - AbAiZpZmZy^ZA = 0. (2.34) 

This makes precise the redundancy or over-completeness of the counterterms 
alluded to above. The same procedure applies to Gy(x), and we conclude that 
the requirement of on-shell 0(a) improvement only determines the combina- 
tions of improvement coefficients b m + b^(ZpZ m )~ 2 , b m — by(ZpZ m ZA)~ 1 Zy, 
and b m — bA(ZpZ Ta Zy) ~ 1 Z&- We emphasize that this redundancy is a generic 
feature of tmQCD, and not linked to special choices for the fields or correla- 
tion functions. In particular we note that the third component of the axial 
variation of any composite field eft has the correct quantum numbers to appear 
as an 0(ap q ) counterterm to (p itself. 

In conclusion, 0(a) improved tmQCD as defined here constitutes a one- 
parameter family of improved theories. In view of practical applications it 
is most convenient to choose b m as the free parameter and set it to some 
numerical value. For reasons to become clear in section 4 our preferred choice 
is b m = — However, in the following we will keep all coefficients as unknowns 
and only make a choice at the very end. In order to define on-shell correlation 
functions which are readily accessible to perturbation theory we will first define 
the Schrodinger functional for tmQCD. It is then straightforward to extend 
the techniques of refs. to tmQCD and study the continuum approach of 
correlation functions derived from the Schrodinger functional. 

3 The Schrodinger functional for tmQCD 

This section follows closely Section 5 of ref. j| and ref. @. The reader will be 
assumed familiar with these references, and we will refer to equations there by 
using the prefix I and II, respectively. 

3.1 Definition of the Schrodinger functional 

To define the Schrodinger functional for twisted mass lattice QCD, it is conve- 
nient to follow refs. PjiOf. The Schrodinger functional is thus obtained as the 
integral kernel of some integer power T/a of the transfer matrix. Its Euclidean 
representation is given by 

Z[p',p',C';p,p,C] = J D[E7]DMD$] e" s ™, (3.1) 
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and is thus considered as a functional of the fields at Euclidean times and T. 
From the structure of the transfer matrix it follows that the boundary condi- 
tions for all fields are the same as in the standard framework. In particular, 
the quark fields satisfy, 

P+^\ X0 =Q = P, P-^\ X0=T = P\ 

(3-2) 

^P-\ xo=0 = p, iiP+\ XQ=T = p, 

with the usual projectors P± = i(l±7o). The gauge field boundary conditions 
are as in eqs.(I.4.1)-(I.4.2) and will not be repeated here. 
The action in eq. (jO|), 

S[U,iP,rP} = S G [U]+S F [U,Tp,iP], (3.3) 

splits into the gauge part (1.4.5) and the quark action, which assumes the same 
form as on the infinite lattice ( |2,2| ). Note that we adopt the same conventions 
as in subsect. 4.2 of ||, in particular the quark and antiquark fields are ex- 
tended to all times by "padding" with zeros, and the covariant derivatives in 
the finite space-time volume now contain the additional phase factors related 
to 6 k , (k = 1,2,3). 



3.2 Renormalization and O(a) improvement 

Renormalizability of the tmQCD Schrodinger functional could be verified 
along the lines of ref. pj. However, this is not necessary as any new coun- 
terterm is expected to be proportional to the twisted mass parameter and is 
therefore at least of mass dimension 4. One therefore expects the Schrodinger 
functional to be finite after renormalization of the mass parameters and the 
gauge coupling as in infinite volume and by scaling the quark and anti- 
quark boundary fields with a common renormalization constant [11|. This 



expectation will be confirmed in the course of the perturbative calculation. 

The structure of the new counterterms at 0(a) is again determined by the 
symmetries. These are the same as in infinite space-time volume, except for 
those which exchange spatial and temporal directions. The improved action, 

S imp r[U, i>] = S[U, fj>] + SSy[U, V>] + SS G>b [U] + SS F>h [U, V], (3-4) 

has the same structure as in the standard framework, in particular, SS V and 
6SQ,b are as given in eqs. (1.5.3) and (1.5.6). The symmetries allow for two 
new fermionic boundary counterterms, 

<D± = ip q ^ 5 T 3 P±^j. (3.5) 
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The equations of motion do not lead to a further reduction and the action 
with the fermionic boundary counterterms at 0(a) is then given by 

5S F>b [U,$,i;] = a 4 ^{(c s -l)[d s (x) + d s '(x)] 

X 

+ (c t -i)[d t (x)-a t / ( x )] 

+ (6 1 -l)[Q 1 (x) + Q' 1 (x)] 
+ (6 2 -l)[Q 2 (x) + Q 2 (x)]}. (3.6) 
Here, we have chosen lattice operators as follows, 

Qi(x) = i^{x) lb T 3 ^{x)\ xo=a , (3.7) 
Qi(x) = i^(x) l5 T 3 i>(x)\ xQ=T _ a , (3.8) 

Q 2 (x) = z/z q p(x)7 5 <r 3 p(x), (3-9) 
Q' 2 (x) = i / u q p , (x) 75 r 3 / 9 , (x), (3.10) 

and the expressions for the lattice operators O s ^ and Og,t are given in eqs. (1.5.21)- 
(1.5.24). Note that the improvement coefficients are the same for both bound- 
aries, as the counterterms are related by a time reflection combined with a 
flavour exchange. 

3.3 Dirac equation and classical solutions 

For Euclidean times < xo < T the lattice Dirac operator and its adjoint are 
formally defined through 

= (D + SD + m + ifi q75 T 3 )^(x), (3.11) 

Sip(x) 

= ^X^ + ^ + ^o + WsA (3.12) 

where SD = 5D V + 5D\> is the sum of the volume and the boundary 0(a) 
counterterms. Eq. (II. 2. 3) for the volume counterterms remains valid, whereas 
for the boundary counterterms one obtains 

SD h ^(x) = (c t -l)i{4 ,a[V'(x)-t/(x-a6,0)- 1 P + ^(x-a6)] 

+ 5 X0 , T -a [i>(x) - U(x, 0)P_ i>{x + a6)] } 

+ (bl ~ 1) [<5 X0 , a + 5 X0 ,T-a] ^q75T 3 V'(x). (3.13) 
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We observe that the net effect of the additional counterterm consists in the 
replacement p q — » b\p q close to the boundaries. Although a boundary 0(a) 
effect is unlikely to have a major impact, we note that the presence of this 
counterterm with a general coefficient b\ invalidates the argument by which 
zero modes of the Wilson-Dirac operator are absent in twisted mass lattice 
QCD. To circumvent this problem we remark that the counterterm may also 
be implemented by explicit insertions into the correlation functions. As every 
insertion comes with a power of o, a single insertion will be sufficient in most 
cases, yielding a result that is equivalent up to terms of 0(a 2 ). 

Given the Dirac operator, the propagator is now defined through 

{D + 5D + m + ifi ql5 T 3 )S(x,y) =a- 4 5 xy , 0<x <T, (3.14) 

and the boundary conditions 

P + S(x,y)\ xo=0 = P„S(x,y)\ X0=T = 0. (3.15) 

Boundary conditions in the second argument follow from the conjugation prop- 
erty, 

S , (x,y) t = 7 5 r 1 S'(y,x)75r 1 , (3.16) 

which is the usual one up to an exchange of the flavour components. 

As in the standard framework [11,3]], it is useful to consider the classical 
solutions of the Dirac equation, 



D + 6D + mo + iw**)M*) = °> ( 3 - 17 ) 
4i(x) fBt + $Bt + m + W5T 3 ) = 0. (3.18) 



Here, the time argument is restricted to < xq < T, while at the boundaries 
the classical solutions are required to satisfy the inhomogeneous boundary 
conditions ( |3.2[) . It is not difficult to obtain the explicit expressions, 

^cl(x) = c t a 3 ^{s(x,y)L/(y-a6,0)- 1 P +/ 9(y)| w=a 
y 

+ S(x,y)U( y ,0)P^p'(y)\ yo=T ^ a }, (3.19) 

rp cl (x) = w 3 J2{p~(y) p -U(y-ab,0)S(y,x)\ yo=a 
y 

+ p'(y)P + U(y,0)- 1 S(y,x)\ yo=T _ a }, (3.20) 

which are again valid for < xq < T. Note that these expressions are exactly 
the same as in ref. Q, except that the quark propagator here is the solution 
of eq. (|3l4|) . 
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3.4 Quark functional integral and basic 2-point functions 



We shall use the same formalism for the quark functional integral as described 
in subsect. II. 2.3. Most of the equations can be taken over literally, in par- 
ticular, eq. (II. 2. 21) holds again. The presence of the twisted mass term 
merely leads to a modification of the improved action of the classical fields, 
[eq. (II. 2. 22)], which is now given by 



,impr 

[U,1p c l,ljj c l] 



J^{&2a/iq [/?(x)i7 5 T 3 /)(x) + p'(x)i-f 5 T 3 p'(x)] 



+ c s a[p(x)7fc^(V fc + Vfc)p(x) 

+ / 9 / (x) 7 4(V fc + V£)p , (x) 



5t [p(x)U (x - aQ, 0)4> c i(x] 
+ p'(x)[/(x,OrVci0r) 



I XQ=a 



x =T- 



-«]}■ (3-21) 



The quark action is a quadratic form in the Grassmann fields, and the func- 
tional integral can be solved explicitly. Therefore, in a fixed gauge field back- 
ground any fermionic correlation function can be expressed in terms of the 
basic two-point functions. Besides the propagator already introduced above, 

[i>(x)i;(y)} F = S(x,y), (3.22) 

we note that the boundary-to-volume correlators can be written in a conve- 
nient way using the classical solutions, 



[C(x)$(j/)] F 

ty(*)C(y)] F 
[C'(xMy)] F 
ty(*)C'(y)] P 



^(x) ' 

Sp(y) ' 



(3.23) 
(3.24) 
(3.25) 
(3.26) 



The explicit expressions in terms of the quark propagator can be easily ob- 
tained from eqs. (3.19) and ( p. 20 ), and coincide with those given in ref. Q 
The boundary-to-boundary correlators can be written as follows, 



[C(x)C'(y)] F 
[C'(x)C(y)] F 



t t P-U(x - a6,0) ty(x)c'(y)] F 
ctP+Uix^)- 1 [V(x)C(y)] F 



xo=a 



xo=T—a 



(3.27) 
(3.28) 
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The correlators of two boundary quark fields at the same boundary receive 
additional contributions due to the new boundary counterterms, viz. 



[C(xR(y)] 



[C'(x)C'(y) 



~2 



P-U(x - a6,0)S(x,y)U(y - 06, O^P+I 



--yo=a 



a~ 2 S 



xy> 



ef P + U(x,0)- 1 S(x,y)U(y,0)P- 



\x =y =T-a 



Pa 



a" 2 <5 



xy 



We finally note that the conjugation property ( 3.16f) implies, 

^(^)C(y)]^ = 75T 1 [C(y)^(x)] F 75T 1 , 

[C(x)C'(y)] F = 75r 1 [C'(y)C(x)] F 7 5 r 1 , 
[C(x)C(y)]i = 7 5 r 1 [C(y)C(x)] F 7 5 r 1 ) 

and analogous equations for the remaining 2-point functions. 
3.5 SF Correlation functions 

With this set-up of the SF we now define a few on-shell correlation functions 
involving the composite fields of Sect. 2. With the boundary source 



(3.29) 

(3.30) 

(3.31) 
(3.32) 
(3.33) 



o a = « 6 Ec"(y)75^ a c(z), 



we define the correlation functions 

rub 



f?(xo) = -(A a (x)O b ), 
f?(x ) = -(P a (x)O b ), 
tf(xo) = -(V a (x)O b ). 



(3.34) 

(3.35) 
(3.36) 
(3.37) 



In the following we restrict the isospin indices to a, b G {1, 2}. It is convenient 

(3.38) 



to define the matrix p2 . IS ] 



H(x) = a 3 Y, 



8ip c \(x) 



s P (y) 



Its hermitian conjugate matrix is given by 
^(x)t = a 3 X)75T 1 



W 75T ' 



(3.39) 
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and the correlation functions can be expressed in terms of H{x), viz. 

fx (*o) = (|tr {H(x)^ 5 T X T l r a H(x)r b r l }) G . (3.40) 

As in ref. Q the bracket (• • • }g means an average over the gauge fields with 
the effective gauge action, 

S cS [U] = S G [U] + 5S Gth [U] - In det (D + SD + m + i/i q 7 5 r 3 ) , (3.41) 

and the trace is over flavour, Dirac and colour indices. The gamma structures 
are Tx = 7o75>75;70) where X stands for A, P and V respectively. 

3.6 Reducing the flavour structure 

In order to carry out the flavour traces we introduce the flavour projectors 

Q± = i(l±r 3 ). (3.42) 
Inserting the flavour decomposition, 

H{x) = H + (x)Q + + F_(x)Q_, (3.43) 



into the expression eq. (3.40) leads to 



/x(*o)= E tr{Q i rVQ i r 6 r 1 }(itr { J ff i (x)t 75 r x ^(x)}) G . (3.44) 

Since we restrict the indices a and b to values in {1,2} this expression further 
simplifies leading to 

fx (*o) = £ tr{Q,rV}(itr {fT i (x)t 76 r x H i (x)}) G . (3.45) 

In order to simplify the expressions further, we now study the behaviour under 
a parity transformation combined with the exchange /i q — ► — /x q . Notice that 
the parity transformation also transforms the background fields, in particular 
it implies 9k — * —0k (k = 1,2,3). On the matrices H±{x) this transformation 
acts according to 

H±(x) — » -yoH^x), (3.46) 

where x = (xq, — x) is the parity transformed space-time argument, and we 
recall that H±{x) depend implicitly on the background gauge field. After 
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averaging over the gauge fields and due to parity invariance of the effective 
gauge action ( 3.41| ) one then finds 



tr [h ± (x)^ 5 T x H ± (x)}) g = r,(X) ( tr {h t (x)^ 5 T x H t (x)}) q , (3.47) 

where the sign factor depends on whether Tx commutes (??(X) = — 1) or anti- 
commutes (r/(X) = 1) with 70. Using this result in eq. (|3.45| ) it follows that 

fl 2 (xo) = f¥(x ) = / v Vo) = 0. (3.48) 

Furthermore, the exact U(l) flavour symmetry implies that 

/x fro) = /x(*o), /x(*o) = -fk 2 (x ), (3.49) 

so that we may restrict attention to the following non-vanishing correlation 
functions: 

/AVo) = 4(tr {h + (x)^ H + (x)}) g , (3.50) 

fi\x ) = ±(tv{H + (x) j <H + (x)}) G , (3.51) 

/v(*o) = |(tr {iT+Cx)^^^)}^. (3.52) 

Note that eq. ( |3.47| ) has allowed to eliminate the dependence upon the second 
flavour component H_(x). This is convenient both for perturbative calcula- 
tions and in the framework of numerical simulations. 

4 O(a) improvement of the free theory 

We determine the improvement coefficients in the free theory, which is obtained 
by setting all gauge links to unity. In this context correlation functions of 
quark and antiquark fields are suitable on-shell quantities which ought to be 
improved. We may therefore consider the improvement of the one-particle 
energies, the quark propagator and basic 2-point functions in the Schrodinger 
functional, in addition to the SF correlation functions introduced in section 3. 



4.1 The free quark propagator 

All correlation functions in the SF are obtainable from the quark propagator, 
which can be computed using standard methods M. We set the standard 
improvement coefficients to their known values [Q], 

c t = c s = 1, (4.1) 
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and compute the propagator assuming b\ = 1. As discussed in sect. 3, any 
other value can be obtained by insertion of the corresponding boundary coun- 
terterm. The propagator can be written in the form 

S(x, y) = (l>t + mo - i/x q75 T 3 ) G(x, y), (4.2) 

where G(x, y) is given by 

G(x,y) = L- 3 Y,^ y) [G + (p,x ,y )P + + G^(p,x ,yo)P-] , (4.3) 



with the functions 

G + (p;x ,y ) = N{p + ){M^{p + ) 



e -u(p+)(\x -y Q \-T) _ e u>(p+)(x +y -T) 



+ M + (p 

G-(p;x ,y ) = G+(p; T - x , T - y ). 



e u(p+)(\x -yo\-T) _ e -u(p+)(x +y -T) 



(4.4) 
(4.5) 



Here, M±(p + ) = M(p + ) ± ipfi (II.3.17), with M(p) as defined in eq. (II.3.6) 
and p^ = p^ + 0^/ L. Furthermore, we recall that in the above formulae it is 
understood that Po = Pq = iui(p + ), where for given spatial momentum q the 
energy w(q) is obtained as the solution of the equation 



sinh 



a 

2«(q) 



a J q 2 + ^1 + (m + |aq 2 
2 1 1 + a(m + iaq 2 ) 



(4.6) 



Finally, using again the notation of ref. 0], the normalization factor is given 
by 



-2ip^A(p + )R{p+)e^ p+ ^ 



(4.7) 



4.2 Improvement conditions and results 



In the free quark theory, the quark energy a; is a suitable on-shell quantity. 
At zero spatial momentum it coincides with the pole mass, which is related to 
the bare masses through 



cosh onip = 1 + ^a 2 (mg + Ma) /(I + amo). 
Up to terms of 0(a 2 ) one then finds (m c = at tree level) 



m, 



+ Mq) (1 - am q ) + 0(a 2 ). 



(4.8) 
(4.9) 
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Replacing the bare masses by the renormalized 0(a) improved mass parame- 
ters and requiring the absence of 0(a) artifacts one obtains 

&m = -i 6 M + 6 m + i = 0, (4.10) 

and the same condition is obtained from the 0(a) improved energy at finite 
spatial momentum. One may wonder whether it is possible to get an additional 
condition by considering the improvement of the quark propagator itself. This 
is not so, for the reasons given in subsection 2.5. As an illustration we consider 
the quark propagator fl4.2|) in the limit of infinite time extent T with the limit 
taken at fixed xq — T/2 and yo — T/2. This eliminates the boundaries both 
at xo = and xq = T, so that one is left with the improvement of the mass 
parameters, and of the quark and antiquark fields, viz. 

ip R = (l + b^arriQ + 6^ia/i q 7 5 r 3 ) ip, (4.11) 

3 



ip R = if) [1 + b^am + b^ia[i q -f 5 T j . (4.12) 

Requiring the quark propagator to be 0(a) improved we find the usual result 
of the untwisted theory, = b^ = |, and 

bj = b^, 2^-6 m -i = 0, 2^ + ^ = 0, (4.13) 

i.e. 3 equations for 4 coefficients. Similarly, by studying the SF correlation 
functions of the improved quark bilinear fields we find the standard results of 
the untwisted theory, ca = cy = and 2b^ = b^ = by = bp = 1, and the 
following conditions involving the new coefficients, 

bi-\(b m + \) = 1, (4.14) 

b„ + ~b m + ± = 0, (4.15) 

bA-(6m + |) = 0, (4.16) 

bv-(b m + \) = 0. (4.17) 

Furthermore, from the O(a) improvement of the basic 2-point functions we 
also obtain 

b 2 = 1. (4.18) 

The fact that b^ and b\ are not determined independently is again due to 
the invariance of the continuum theory under axial rotations of the fields and 
a compensating change in the mass parameters. Hence our findings in the 
free theory are completely in line with the general expectation expressed in 
subsect. 2.5. Choosing b m as the free parameter and setting it to — \ leads 
to 6 M = 6a = by = and b\ = 1, while e.g. for 6 m = the tree level value 
b\ = 5/4 is somewhat inconvenient. 
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5 The one-loop computation 

We now want to expand the correlation functions to one-loop order. We work 
with vanishing boundary values C k and Ci. The gauge fixing procedure then 
is the same as in ref. Q and will not be described here. In the following we 
only describe those aspects that are new and otherwise assume the reader to 
be familiar with refs. [§|J5]]. 



5.1 Renormalized amplitudes 

Once the flavour traces have been taken, the one-loop calculation at fixed 
lattice size is almost identical to the standard case @||. In order to take the 
continuum limit at fixed physical space-time volume, we then keep m R , fi K , xq 
and T fixed in units of L. Here the renormalized mass parameters are defined 
in a mass-independent renormalization scheme which may remain unspecified 
for the moment. 

To first order of perturbation theory the substitutions for the coupling 
constant and the quark mass then amount to 



9l 
m 



4 ) 



al + o{ g 



(5.1) 
(5.2) 
(5.3) 



where the precise form of the coefficients 



■in, 



(0) 



1 — yT- 2am R — a 2 /i R 



(5.4) 



m, 



(i) 



(i) 



in 



+ a/4 



z (i) 



6«< 



a m 



(0) 







(0) 



am. 



(0) 




(5.5) 



u (0) 



Mr j 



-u(°) 

Mq 



{4" 



+ &!,'W< 0) } 



(5.6) 
(5.7) 



is a direct consequence of the definitions made in subsect. 2.1, and already 
includes the tree-level results obtained in the preceding section with the par- 



ticular choice b 



(0) 



l 

2" 
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The renormalized correlation functions, 

[fy(xo)]n = Zy(l + byam q )Z 2 (l + 6(am q ) 2 

x {/^(x ) + M/Va(*o)} , (5.8) 
[f^(x )] R = Zp(l + bpam q )Z^(l + b ( am q ) 2 f^(xo), (5.9) 
[/a^o^r = ^a(1 + b A am q )Z 2 (l + 6 c am q ) 2 

x {/a(*o) + CAaSo/pVo) - 6 A a/i q /v fco)} , (5-10) 

have a well-defined perturbation expansion in the renormalized coupling g R , 
with coefficients that are computable functions of a/L. For instance the ex- 
pansion of [/v 2 ]r reads 

[fy 2 (x )] R = /?(xo) (0) + 5 l{/?(xo) (1) +m«^-/^(xo)(°) 
+ (4 1} + 2Z« + am R [6« + 26^]) /^ 2 (^o) (0) 
+ ^ 1) ^/v(*o) (0) +a^bMfl\x )®}, (5.11) 

where terms of order a 2 and g R have been neglected, and it is understood that 
the correlation functions are evaluated at mo = and fi q = [x q ^ . 

Following ref. we now set xq = T/2 and scale all dimensionful quantities 
in units of L. With the parameters z m = m R L, z^ = fi R L and r = T/L we 
then consider the dimensionless functions, 



h A (0,z m , 




a/L) 


= [/a hjk I0=t/2 , 


(5.12) 


hy(6, z m , 


Zfi,T, 


a/L) 


= [/v 2 ( a; o)]R| a . 0=T / 2) 


(5.13) 


hp (9, Z m , 


Zfi, T, 


a/L) 


= [fp( x o)}n L =T/2' 


(5.14) 


hdA(9, Zm, 


ZfM,T, 


a/L) 


= 4[/A^)k IO= T/2' 


(5.15) 


hay (6, Z m , 


Zfn T, 


a/L) 


= Ld [fy 2 (x )] R \ XQ=T/r 


(5.16) 
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One then infers, 



+ (zW + 2Z« + |z m [6W + 26[ 1 )]), ^ 



Zip = 



/idV — 



u W = 



= - Z JL 



L L 2 

A 4 A* £ 



(5.17) 



+ (4 1 ) + 2Z« + ^ m [6« + 26«]) go }, (5.18) 

«o + + cf )n 2 + am^u 3 + anvils + z^ 1 ^ 

+ (Z P X) + 2Z« + + 26«]) Uo }, (5.19) 

hdA = wo + 5||wi + c[ 1 ^ 2 + am^'u;3 + c^tt;4 + a j uWi(;5 



(5.20) 



+ (4 1} + 2Z« + ^ m [6« + 26«]) r }. (5.21) 
Since we are neglecting terms of order a 2 , the expansions, 



(5.22) 
(5.23) 



_L L ±y J 

may be inserted in Eqs. (5.1 <" )-( 5.21| ). All the coefficients Vi, . . . ,rj are still 
functions of t ? i i z m cind 2^. Analytic expressions Ccin be derived for those 
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coefficients involving the tree level correlation functions or the 0(a) countert- 
erms. Their asymptotic expansions for a/L — > are collected in Appendix B. 
The coefficients vx,...,r± are only obtained numerically and definite choices 
for the parameters had to be made. We generated numerical data for 9 = 
and 9 = 0.5 for both T = L and T = 2L and various combinations of the 
mass parameters z m and ^ with values between and 1.5. With these 
parameter choices the Feynman diagrams were then evaluated numerically in 
64 bit precision arithmetic for a sequence of lattice sizes ranging from L/a = 4 
to L/a = 32 (and in some cases to L/a = 36). 



5.2 Analysis and results 

The renormalization constants are determined by requiring the renormalized 
amplitudes to be finite in the continuum limit, and by the requirement that 
the tmQCD Ward identities be satisfied [0J. A linear divergence is cancelled 
in all amplitudes by inserting the usual one-loop coefficient arni \ or equiv- 
alently a series which converges to this coefficient in the limit a/L — > Q. 
We choose the lattice minimal-subtraction scheme to renormalize the pseudo- 
scalar density and the quark boundary fields, and the one-loop coefficients are 
then given by [with C F = (N 2 - l)/2N], 

Z^ = -^HL/a), 2zf = -Zf. (5.24) 

The current renormalization constants, and the renormalization of the stan- 
dard and twisted mass parameters are determined by the Ward identities. For 
the one- loop coefficients we expect [14,15,0], 



= -0.087344(2) C F , (5.25) 

Z$p = -0.097072(2) C F , (5.26) 

Z^ = -Z p 1} - 0.019458(1) C P , (5.27) 

ZjP = "4 1} - (5-28) 

With our data we were able to compute the one-loop coefficients of the combi- 
nations Z m Zp/ZA and Z^Zp/Zy, as well as the logarithmically divergent parts 
of all one-loop coefficients. Complete consistency with the above expectations 
was found, and we shall adopt these results in the following. 

The corresponding coefficients in other schemes differ from those above 
by a- independent terms. With the renormalization constants chosen in this 
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way we find e.g. for the combination of separately diverging terms appearing 
in the curly bracket of fl5.19| ) 



ui + am^us + 



(4 1) +2ZJ 1) )uo 



(-1) 

m "3 



7(1), „ 



(-1) 



U +U 1T + O(a 2 /L 2 ) 



(5.29) 



where ZYj are functions of r,9,z m and and ^ are coefficients of L/a 
in the expansion of itj for L/a — > oo. Evidently similar equations hold for 
the other functions ui, gi, it?i, n. It is important to note that we expect no 
terms involving (a/L) ln(L/a) on the right hand side of ( |5.29| ) because we 
have imposed tree level improvement, and this was indeed seen in our data 



analysis. Moreover there are no terms ~ Z^a/L or ~ Z^'a/L on the left 
hand side above because of Eq. ( B,23j ); thus the coefficient U\ is (contrary to 
IAq) independent of the renormalization scheme. Estimates for the coefficients 
Ui,Vi,--. were obtained for the various data sequences using the methods 
described in [16]. 

Now the improvement coefficients are determined by demanding that the 
renormalized amplitudes approach the continuum limit with corrections of 
0(a 2 /L 2 ). For the cancellation of the 0(a) terms the following equations 
should be satisfied (for undefined notation see Appendix B): 



b {1) u {1) 

°1 "6 



z ^ u ia '3 



(\(») 
V ^0 



°1 '6 



Vi+Vi, (5.30) 
Qi + Qi, (5.31) 
Wi+Wi, (5.32) 
Wi + Wi,(5.33) 
IZi + K 1 . (5.34) 



In these equations all terms involving improvement coefficients which are nec- 
essary also in the untwisted theory, have been collected in the terms U\, . . . on 
the right hand sides and they are specified in equations ( p,45| )-( |B.49 ). The nu- 
merical values of these improvement coefficients, obtained in previous analyses 
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!§, are: 



c? = 


-0.01346(1) C F , 


(5.35) 


c ( P = 


—0 005680(2) C-p 


(5.36) 


»? = 


-0.06738(4) C F , 


(5.37) 


"m 


-0.07217(2) C F , 


(5.38) 


°A — 


0.11414(4) C F , 


(5.39) 


6 {1) - 

Oy — 


0.11492(4) C F , 


(5.40) 


Op — 


0.11484(4) C F . 


(5.41) 



Before we proceed with the numerical analysis of equations fl5.30| )-( 5.34 ), 
it is essential to note that using the identities ( |B.38|) -(|B.42i) they can be rewrit- 
ten as 



z //(V~ 1} + 6 /{1) o (0) 
z m o^ v 5 -f- o A q 



y {1) v {1) 



V V 



= Vi+Vi, 
= Qi + Qi, 
= Wi+Si, 
= Wi + Wi, 
= Tli + Hi, 

where the primed coefficients appearing here are defined through 



z '(V _1) 

z m> ^5 + °A r °1 w 6 
til) J- 1 ) 



Oy Wq 



°1 '6 



(5.42) 
(5.43) 
(5.44) 
(5.45) 
(5.46) 

(5.47) 
(5.48) 
(5.49) 
(5.50) 



In other words, from our equations we can only obtain information on four 
linearly independent combinations of the new improvement coefficients ap- 
pearing in the twisted theory. This was in fact to be anticipated from our 
general discussion in subsect. 2.5, where we argued that we are free to chose 
for example the coefficient b$ as we please. 



bp 


= &?> 


+ 6« 


bp 


= <>? 


- ISM 


yW 

°A 


- b {1) 

— °A 




y(i) 
°v 


- b W 
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Since our equations are over-determined and also having generated such 
a large selection of data sets, we had many ways to proceed to determine 
the coefficients b'fi , t>P , b~P and b~P , and a multitude of consistency checks 
on the results. We first note that if we consider the linear combination of 
amplitudes h^A — 2z m /ip and hay + 2z^hp associated with the PCAC and 
PCVC relations respectively we obtain 

- 2zluf ) ~bf ) = Wi + Wi - 2z m {Ui + Ux) , (5.51) 
-2 Zli z m uf (bp + bp) = 1Z 1 +H 1 + 2z^ (Ui + U x ) . (5.52) 



With knowledge of the right hand sides, each equation determines a partic- 
ular linear combination of improvement coefficients. In these equations the 
boundary coefficient b\ does not appear as expected. On the other hand the 



coefficient b'P is all that appears on the left hand sides of Eqs. (5.44),(5.4£ 
for the data sets with z m = 0. 



By solving simultaneously the three equations fl5.42| ), ( |5.44[ ) and (5.45) for 



one data set with z m ^ 0, we could obtain the three coefficients^ b'P , b'P and 
b'P (and of course analogously for the equations involving the vector current). 



We also extracted the two coefficients b'ji by solving just Eq. ( 5.44[ ) for 

two different data sets (of which at least one has z m ^ 0). 

Unfortunately due to rounding errors, the one-loop cutoff effects like U\ 
were rarely determined better than to within a few percent. The consequence 
of this was that many routes of analyses described above and when applied 
to various (combinations of) data sets, led to results for the improvement 
coefficients with very large errors. Nevertheless there remained sufficiently 
many analyses which delivered useful results with relatively small errors, and 
in these cases all results were consistent with each other and with our following 
"best estimates": 

bp = -0.103(3) C F , (5.53) 
bp = 0.035(2) C F , (5.54) 
bp = 0.086(4) C F , (5.55) 
bp = 0.074(3) Op. (5.56) 

As one practical choice for applications in numerical simulations we ad- 
vocate b m = — | to all orders of perturbation theory, which would result in 
setting b~m = in the above equations. 

Particularly good results were obtained e.g. with the data set z m = , z M = 0.5 ,8 — 0, 
where we in fact had data up to L/a = 36. 



23 



6 Conclusions 



In this paper we have introduced the set-up of 0(a) improved twisted mass 
lattice QCD in its simplest form with two mass-degenerate quarks. In per- 
turbation theory to one- loop order we have verified that 0(a) improvement 
works out as expected. We have identified the new counterterms and com- 
puted their coefficients at the tree- level and to one- loop order. In practice 
perturbative estimates may be satisfactory, as tmQCD has been primarily 
designed to explore the chiral region of QCD, where the contribution of the 
new counterterms should be small anyway. This expectation is confirmed by 
a non-perturbative scaling test in a physically small volume, which employs 
the perturbative values of the new improvement coefficients reported here [19|. 
However, a non-perturbative determination of some of the new coefficients is 
certainly desirable and may be possible along the lines of ref. ||. 

An interesting aspect of 0(a) improved tmQCD is the absence of any 
new counterterm corresponding to a rescaling of the bare coupling go- This 
singles out the choice for the angle a = ir/2 for which the physical quark mass 
is entirely defined in terms of the twisted mass parameter. A quark mass 
dependent rescaling of go is hence completely avoided, and one may hope 
that this eases the chiral extrapolation or interpolation of numerical simula- 
tion data. Furthermore, using the over-completeness of the counterterms (cf. 
subsect. 2.5) to fix b m exactly, no tuning is necessary to obtain a = ir/2 up 
to 0(a 2 ) effects, provided the standard critical mass m c and the standard 
improvement coefficients of the massless theory c sw and ca are known. We 
also note that, at a = tt/2, both sides of the exact PCVC relation are au- 
tomatically renormalized and 0(a) improved. This can be exploited for an 
0(a) improved determination of F n |20|| , as the vector current at a = tt/2 is 
physically interpreted as the axial current 0. 

In the future one may wish to extend the framework of 0(a) improved 
tmQCD to include the heavier quarks in the way suggested in ref. [Q. The 
analysis of 0(a) counterterms still remains to be done, but we do not expect 
any new conceptual problems here. 

Finally, we have defined the Schrodinger functional for tmQCD, based 
on the appropriate generalisation of Liischer's transfer matrix construction 
for tmQCD. We expect that the Schrodinger functional will be useful in the 
determination of hadronic matrix elements along the lines of refs. [17,18], and 
work in this direction is currently in progress [pQ , pi]1 . 

This work is part of the ALPHA collaboration research programme. We 
are grateful to P.A. Grassi for discussions and his collaboration in the tmQCD 
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project. Thanks also go to M. Liischer, R. Sommer and A. Vladikas for use- 
ful comments and discussions. S. Sint acknowledges partial support by the 
European Commission under grant No. FMBICT972442. 



A The transfer matrix for twisted mass lattice QCD 

In this appendix we briefly indicate the generalization of the transfer matrix 
construction for twisted mass lattice QCD with c sw = 0. We use the original 
notation of ref. || with the conventions of ref. fpLOtl . The transfer matrix as 
an operator in Fock space and as an integral kernel with respect to the gauge 
fields has the structure 

T [U, U'] = fi(U)K [U, U']f F (U% (A.l) 

with pure gauge kernel Kq and the fermionic part 

f F (U) = det(2 K i?) 1 / 4 eMx j P-Cx) exp(-x t 7oMx)- (A.2) 

Here, the operators Xi(x) are canonical (i is a shorthand for colour, spin and 
flavour indices) viz. 

{ Xi (x), x J(y)} = ^«- 3 5 xy , (A.3) 

and B and C are matrix representations of the difference operators 

3 



B = 1 -6K-a 2 K^VfcV fe , (A.4) 

k=l 

3 

C = a^74(V fc + Vt)+m^ q75 r 3 . (A.5) 

k=l 

As in the standard case the positivity of the transfer matrix hinges on the 
positivity of the matrix B, which is guaranteed for < 1/6. This is the 
standard bound which also ensures that the matrix M, 

M = \\n{\Bn- 1 ) , (A.6) 

is well-defined. No restriction applies to the twisted mass parameter, except 
that /z q must be real for the transfer matrix ( A.l ) to reproduce the twisted 
mass lattice QCD action. 



25 



B Analytic expressions for expansion coefficients 

In this appendix we provide explicit analytic expressions for the tree-level 
amplitudes and the counterterms appearing in eqs. ( 5.17] ) — ( |5.21] ) which are 
needed to compute the one-loop amplitudes up to terms of 0(a 2 ). We have 
checked that the analytic expressions correctly reproduce the numerical values 
obtained by directly programming the correlation functions and counterterm 
insertions. 

First we define 



u = ^zl + W^ + zl, (B.l) 

co = cosh(a;r) , (B.2) 

si = sinh(a;r) , (B.3) 

p = uco + 2: m si, (B-4) 

v = Lusi + z m co , (B.5) 

where r = T/L. Then we have uq = + 0(a 2 /L 2 ) etc. with 



Nuj 

3 

P 

32 i „2 



'o 

(0) N(36' 2 + z* + z m v) 



P 2 



For the boundary terms we define 



and then U2 = avlp /L + 0(a 2 /L 2 ) etc. with 



(B.6) 
(B.7) 



(0) Nz^-Zm + V) , > 

% ~ ^2 ' 

4 0) = 2z m uf , (B.9) 

rj 0) = -2^4 0) . (B.10) 



2 3# 2 + z 2 si 

z m + , (B.ll 

P 



U W = 2fuf , (B.12) 

„(., _ 2 ^,_4^ + |)(-„„ + ,) | (B13) 

,<» = 2r,y" - i I^a^± 4 + f^j , (B . 14) 

4 1 ' = 2fw ( ®\ (B.15) 

r W = 2frf ] . (B.16) 
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Similarly, uq = aujp /L + 0(a 2 /L 2 ) etc. with 
„(i) _ 2z M si 7/ (o) 
(i) 2z M si (o) 2Nuz^(-z m + z/) 

(i) 2,^1 (o) 2jVoj (^ + z 2 (l - co)) 

% - —% + ~ 3 - 

«,« = -^4 0) , 
_ 2 ^ si r (o) 

For the derivatives with respect to the mass parameters we have, 
Ui = {L/a)u { r 1] + uf ] + 0(a/L), (i = 3, 5) 

with 

(0) _ (-i) (0) _ (-1) 

and analogous equations hold in all other cases. Defining 

v(l + z m r) 



X 
Y 
X 
Y 



UJ 

p{\ + z m r) 

UJ 

ZjAyr_ + co) 

UJ 

z fl (pr+ si) 



one has 



(-1) , 
u 3 = — + > 

(-1) 2Xvf ] N(v + z m Y) 

(0) 



?3 



2Xg^ JV^(l-y) 



P P 2 



r M) - -2z u M) 
'3 ^z fl u 3 , 
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and 



( 

4" 




_ yl "0 

p 


4- 


-1) 


2Xv® 


4- 


-1) 


2X$ 
P 


w 5 


-1) 


- 2z u { ~ 1] 


r (- 
r 5 


-1) 


= -2^4" 



+ 



+ 



ujp 

N{2z^ + z m Y) 
N{z^Y -z m + v) 



2u, 



(o) 



Note the identities 



= 


2z^u^ 


-i) 


-2z u { - 1} 


-«?>, 


= 


2Z(j,vf 


1) 


-2z v { - 1] 


-4 1} + ^\ 


= 


2z li qf 


i) 


— 2z m qi i 




= 


2z /Jl w$ 


-l) 


— 2z m wi i 


" 4 1} + 2rf 


= 


2z il r\ > 


i) 


— ?7 r^ 1 ^ 


- r« - 24 0) 



(B.33) 
(B.34) 

(B.35) 

(B.36) 
(B.37) 

(B.38) 
(B.39) 
(B.40) 
(B.41) 
(B.42) 



The remaining coefficients to be specified are V4 = av^ /L + 0(a 2 /L 2 ) and 
w 4 = aw^/L + 0(a 2 /L 2 ) with 



(1) 



2Nuj 2 u 
4Nlo 3 



(B.43) 
(B.44) 



Finally we specify the terms IA\ , . . . appearing on the right hand side of 



(B.45) 
(B.46) 
(B.47) 



equations (5.30)^(5.34): 




Vi 




-zlbg^ + z m [b^ + 2bf}v^+c^v 4 , 


Si 


= 


-zlb^ q ( f l) +z m [b^+2bf] q ^, 


Ui 


= « } 


-zl l b^u^+z m [b^+2bf]uf, 


m 


= ciM } 


- z 2 m bg) W ^ + z m [bW + 26«]4 0) + ofw, 






— 2^j6^V 3 ^ + z m [6y^ + 26^ ^]tq ^ . 



(B.49) 
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